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出 版 前 言 


2011 年 12 月 11 日 是 西安 交通 大 学 杰出 校友 钱学森 先生 
的 百年 诞 展 。 为 缅怀 钱学森 学 长 ,学 习 他 的 科学 思想 和 卓越 风 
范 ,展示 其 丰功伟绩 和 人 格 魅 力 , 西 安 交 通 大 学 举办 了 “纪念 钱 
αξίες 100 周年 "系列 活动 :作为 制 片 方 之 一 ,参与 西部 电影 
集团 摄制 传记 故事 片 (钱学森 ): 与 中 央 电 视 台 合作 ,出 品 纪录 片 
《实验 班 的 故事 一 一 沿 着 钱学森 走 过 的 路 ); 扩 建 钱学森 生平 业 
绩 展 馆 .向 校内 外 开放 ;举办 钱学森 科学 与 教育 思想 研讨 会 ;出 
版 发 行 (钱学森 力学 手稿 )《 钱 学 森 年 谱 ( 初 编 )》《 钱 学 森 第 六 
次 产业 革命 思想 探 微 丛 书 ) 等 。 

钱学森 先生 在 美国 深造 和 工作 期 间 留 下 大 量 珍贵 手稿 ,这 
些 手稿 真实 展示 了 钱学森 先生 博大 精深 的 学 识 . 开 拓 求 实 的 精 
神 和 严谨 奋进 的 作风 ,是 钱 老 勇 攀 科 学 高 峰 和 严谨 治学 的 集中 
体现 ,这 里 ,我 们 将 部 分 原稿 整理 汇集 成 册 . 出 版 (钱学森 力学 
手稿 ,作为 钱 老 百 年 诞 展 的 献礼 。 

人 钱学森 力学 手稿 ) 共 10 卷 .包含 两 部 分 内 容 。 第 一 部 分 是 
草稿 .包括 扁 壳 、 球 过 和 圆柱 这 屈曲 分 析 的 公式 推导 和 数值 演 
Я. 在 研究 圆柱 壳 轴 压 届 曲 问题 时 ,为 了 求 得 圆柱 壳 体 的 临界 
压力 .在 有 关 的 五 百 多 页 草稿 中 ,对 多 达 二 十 多 种 可 能 的 屈曲 模 


态 逐 一 进行 公式 推演 和 数值 计算 ,最 终 才 找到 满意 的 并 在 论文 
中 采用 的 屈曲 模 态 。 仔 细 观 察 草 稿 中 的 数据 列表 ,每 个 数字 有 
效 位 数 都 长 达 八 位 ,在 手 摇 机 械 式 计算 机 作为 主要 计算 工具 的 
年 代 , 这 串 串 数字 凝聚 着 多 少 现今 难以 想象 的 艰辛 劳动 。 

第 二 部 分 是 手稿 ,以 航空 航天 工程 为 核心 ,涵盖 空气 动力 
学 、 固 体力 学、 火箭 技术 ,工程 控制 论 和 物理 力学 等 领域 的 部 分 | 
学 术 论 文 手稿 ,打印 稿 和 讲义 。 

《钱学森 力学 手稿 ) 是 在 西安 交通 大 学 校 领导 的 大 力 支持 
下 ,由 西安 交通 大 学 航天 航空 学 院 沈 亚 鹏 教授 整理 完成 。 图 书 
出 版 过 程 中 得 到 了 西安 交通 大 学 党 委 宣 传 部 .校友 关系 发 展 部 、 
图 书馆 .航天 航空 学 院 等 的 积极 协助 ,在 此 深 表 感谢 。 
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Heue-shen Tsien 
Guggenheim Aeroanutical Laboratory 
Massachusetts Institute of Technology 


Wind-tunnels are, perhaps, the most useful tool in aerodynamic in- 
vestigations and certainly have contributed much in the modern develop- 


ment of fluid mechanics, It is thus natural, when one turns to a new field 
ef aerodynamics, the aerodynamics of rarefied gases or вирегавгодупаві св, 
that one should think of using the wind tunnel again, Only here it has to 
be adopted to entirely new circumstances and many new problems, both in its 
design and in its operation, appear. It is the purpose of this paper to dis- 
cuss some of these problems, so as to gain an orientation in the new feld 


of experimentation. I 


1. Tunnel Design 

To test models in the wind tunnel at ita test section, it 1s of 
primary importance to obtain a uniform stream at the desired temperature, 
pressure and velocity. For subsonio wind tunnels with ordinary pressures, 
this oan be achieved without much difficulty. For supersonic wind tunnels 
at ordinary pressure, the expansion pert of the tunnel before the test sec- 
tion or the nossle ie first designed to obtain а uniform stream at its exit 
without considering the effects of the viscosity of air. Then the boundary 
layer along the well of the nozzle is calculated with the pressure gradient 
thus determined, Finally the displacement thickness of the boundary layer, 
or the needed space for the boundary layer flow at lower velocity is provid- 
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oaloulated amount. (Ref. 1) This design procedure is found to give satis- 
factory nossles for supersonic wind tunnels. 

However when one tries to use the same design procedure to the super- 
aerodynamics wing tunnel, one is immediately confronted with the difficulty 
of extremely large viscous effect, In other words, the boundary layer will 
be so thick as to occupy the main portion of the nossle passage. То demon- 
strate this effect, let the length of the test section be L and tne width of 
the square test section be b. Then the Reynolds number based upon the oondi- 


tions in the test section is Re = И where U is the velocity in the test seo- 


tiom, If, ав а rough estimate, we take the thickness of the boundary layer 
to be zero at the beginning of the test section and equal to в value calcu- 
lated by the well-known Blasius formula for а flat plate at the end of the 
test section, then 

$ = 2.65 /. Q) 


= 5. %. (2) 


On the other hand, the ratio of the mean free path | and the boundary layer 
thickness Ü is known (Ref, 2) to be equal to 


L М. (3) 
{- un 
where ) is the ratio of specific heats and can be taken as 1.4, М is the Mach 
number іп the test section. By combining (2) and (3), we have 


H Τὴ = оға 


CE (GM E к dee 


and L/b = 2, the boundary layer will completely fill up the test section, 
if the mean free path is equal to 5.6% of the boundary layer thiokness or 
2.8% of the tunnel width, This means that the extremely strong viscous ef- 
fect at low derivities makes the ordinary concept of designing а wind tunnel 
totally inapplicable. 
The extremely thick boundary layer where the velocity increases from а 
small value near the wall to soms supersonic velocity at the center of the 
mossle, also gives subsonio velocities in a rather large portion of the nos- 
Since pressure disturbances downstream can be transmitted upstream in 
subsonic flows, the flow in the test section of ^ low pressure tunnel will 


ale. 


be sensitive to changes in the diffuse even if the main stream velocity at 
the center of the noasle is supersonic. This is of course а new phenomenon 
is superacrodynamics tunnel not found in conventional supersonic wind tunnels, 
The large viscous effects can be also demonstrated by calculating the 
Tatio of frictional loss on the walls of the test section and the shock loss 
in the diffuser after the test section, Consider the diffuser to be a straight 
tube of approximately the same cross-sectional area as the test seotion, then 


ats 


Taking ( to be Hlasive value or ("ΠΣΕ we have 


1.328 (5) 


σα 


ΑΛ zeu t) 


Now if we Se QUIE then the pressure by 


ideal isentropic compression in the diffuser is 
κ fae өнер керекке ы 
mal shock without further recovery, then the actual pressure rise is 


Therefore the pressure loss due to shock 


By combining (5) and (6), the ratio of these two pressure losses {8 


Introducing the moun free path ratio given by (3), жо tave 


(8) 


This relation is plotted in Fig. 2. Therefore if the Mach number ів 2, 
and / = 2 as before, then when the ratio () is 0.056, the ratio of - 
frictional loss to shook loss is 0,628. Hence the frictional loss and the V 
shook loss is of the same order of magnitude, 

These large viscous effects are fully confirmed ty the recent teste on 
the 1" x 1" low pressure wind tunnel at the University of California," The 
test nosale (Fig. 3) was designed for Mach number without considering 
the viscous effect of the medium. During test, the static pressure on the 


* Experimental work done under contact with the Office of Naval Research. 
The author is deeply indebted to Professors F. G, Folsom and E. D. Kane 
for permission to use their unpublished results. 


wall at the exit 


microns" for the two testa presented in Figs. 4 and 5. The apparant Mach 
number "И" is the Mach number oaloulated from the dynamic pressure measured 
by а Pitot tube by using the Rayleigh formula, Since there is the compli- 
oatiom of large viscous effect in the Pitot tube reading as shown in the fol- 
lowing section, this apparant Mach number is only quantitative and cannot be 
taken as the exact value. However it is apparant from Figs. 4 and 5 that the 
boundary layer in the test section 1s indeed very thick, and filla up the 
whole spice, This large boundary layor thickness mikes this spac available 
for th: expansion of the central potential flow if 1t existe, very small. 
Therefore tho maximum Mach number reached at the center of the nomale 4 very 
much smaller than the design Mach number of 4. At the lower prossure, the in- 
fluence of slip nt tho wall is also evident, This hns the tendency to make 
the flow more uniform, However the very low Mach number at the test section 
indioates again the strong viscous effect in converting much of the pressure 
energy into heat energy. 

These elementary calculations and preliminary test results makes it clear 
thet for the design of the nossle and test section for s superaerodynamios wind 
tunnel, it is no longer possitle to separate the compressibility effects and 
the viscous effect. In fact, the concept of boundary layer is also of doubtful 
value due to the extremely small Reynolds number encountered, Therefore to aot- 
ually design such а nossle to obtain the nearest approximation to the ideal uni- 
form flow it will be necessary to use the exact Navier-Stokes equations instead 
of the approximate boundary layer equations. Of course, it may be argued that 
for superaerodynamics, the Navier-Stokes equations for no more exact and addi- 
tional corrections must be added (Ref, 2). However, recent investigations by 


6 
* 1000 microns = 1 mm Hg one atmosphere = 0,760 x 10 "torons. 
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Ап case of slip-flows concerned here and will not essentially alter the flow 
pattern. Henoe for & first approximation just like the non-visoous isentro- 
pic flow as а first approximation for ordinary supersonic nossles, we сап use 
the Mavier-Stokes equations, The simplest case to be considered is certainly 
the axially symmetric nossles. If x is the coordinate in the axial direction, 

the coordinate in the radial direction and u and v are the corresponding 


velocity components, the equations are: 


ж dessipation function 


= stresses tensor 


(9), (20), (12) and (12) together with the equation of states 


-.. (13) 
then determines the five unknown u, v, , and Т. Of course, the actual process 
of making this calculation will be extremely tedious and some approximation 
method of solution may have to be developed. Опе possibility would be to adopt 


the differential equations once with respect to г and thus only try to sat- 


isfy the equations "on the average" over the cross-section of the nozzle. 
The "distribution" of u, v over the cross-section will then be set in the 
form of а polynomic.l in г. Initial study in this procedure is already made 
by B. А. Schaaf (Ref. 4) at the suggestion of the author. 

For ordinary supersonic diffuser, high efficiency of pressure recovery 
Gan be generally achieved by using a long diffuser. However, for superaero- 
dynamios wind tunnel, due to the extremely large loss through friction, long 
diffusers are undesirable. In fact, tho pressure loss сап be reduced by us- 
ing а shortest possible diffuser, 


2. Flow Measurement 
The quantities which determine the fjow field are two out of 
the three variables , , T and the velocity components. The quantities 
„ » Tare related by the equation of states and therefore only two ia nec- 
essary for the determination of all three. Generally for wind tunnel work, 
the quantities actually measured are , „ and , the magnitude of the velo- 
city. 

For the measurement of pressure, а nanometer is used. For ordinary pres- 
sure, one uses а fluid manometer filled with water, alcohol or meroury, How- 
ever for the extremely low pressure encountered in the superaerodynamic flow, 
some other form of manometer is necessary. One of the most successful type is 
the Pirani gage. Тһе conventional form of Pirani gages has n pressure sensiti- 
жібу of about 10? micron. It utilises the change of temperature of a wire 
heated with constant energy caused by а change in the pressure of the gas sur- 
rounding it. The temperature change is measured by the change in the resistance 


of the wire, Тһе wire is looated in a small chamber which is connected to the 


5 


із to be measured, The question of best design of the connecting tube for 
quick response is studies by S. A. Schaaf. (Ref. 5). 


To measure the density , the conventional method utilises the differ- 


ence in the velocity of the light rays in mediums of different density. With 


different optical arrangements, we have the shadowgraph method, the schlieren 
method and the interferometer method. However, if the density of the medium 
is very low as the case of superaerodynamic flow, the sensitivity of these 
methods become extremely poor. For instance in case of schlieren method, the 
percentage change in illumination I by passing through a region of thickness 
b is given by 


Q4) 


where is the air density at 32°F and 1 atmosphere pressure, and 

is the density gradient normal to the light ray. end аге the local length 
and the normal unobscured width of the light source image perpendicular to the 
knife edge. Ж is factor of order 1, determined by the particular optical path 
used, Therefore the sensitivity of the schlieren method decreases with the fac- 
tor ( ). Some improvement сап be made by altering the quantities f and , 
but prmciical limitations and diffraction difficulties do not allow the increase 
of sensitivity to satisfactory values, 

А new approach to this problem density measurement is the method of absorp- Ν 
tion. It is found for instance, that oxygen at low pressures shows a strong ab- 
t sorption band at wave lengths around 1470? А or ultra-violet light. The per- 
| centage absorption is proportional to the number of molecules that meet the light 

тау and is, therefore, proportional to the density of the gas. The measurement 
is then similar to that of the interferometer method where the density is deter- 
mined, A similar method is the utilisation of the after-flow of nitrogen. These 


„„ 


The conventional method for the measurement of velocity is through the 
use of dynamic pressure rise in ^ Pitot-tube, A straight forward application 
of this method is, however, difficult for rarefied gases, The formula used 


is, however, based upon the neglection of viscosity effects. But for rarefied 


gases the viscosity effect ів ef great importance as pointed out in the previ- 
ous section, Then the dynamic pressure would be quite different than that 
given by the usual formula, To estimate this effect, lot us consider the case 
of low Mach number so that compressibility effects can be neglected. Then as 
a firet approximation, take the flow field around the Pitot-tube as that of а 
source of strength in non-viscous flow of uniform velocity 0. The "radius" 
of the tube & is 
.. 


and the stagnation point is located at 


Тһе velocity introduced by the source is then 


By calculating tho viscous stress from this approximate disturbance velocity, 
wo havo for flow slong the axis 


(16) 


Hence ff in the stagnation pressure and 


ο 


number of the Pitot-tube could be of the order of unity. Then the dynamic 

pressure rise ~ is not the usual values $ but а value much lees than 
h that, In faot previous investigation by Baker (Ref. 7) and F. Homann (Ref. 8) 
indicate that the Reynolds number 


where а! is the radius of the mouth of. 
tube, must exceed 30 in order to reach the usual dynamic pressure rise ф 
When the velocity of flow is high, we have the added complication due to 
the shook. Тһе conventional Rayleigh formula for Pitot tubes is supersonic 
flow is based upon the assumption of very thin shock wave ahead of the Pitot 
| tube, Now the thickness of the shook is proportional to the mean free path. 
1 Hence in rarefied flows, the thickness of the shook will be so increased as to 

cause interference with flow in the neighborhood of the Pitot-tube. This to- 

gether with the viscous eífect mentioned in previous paragraph definitely show 

the inapplicability of the Rayleigh formula for supersonic velocity of rarefied 


With the great complications in applying the conventional velo- 

} city measuring device to superaerodynamic flows, опе is naturally led to the 
thought of other avenues of approach. One possibility is the use of hot-wire, 

If the wire diameter is of the order of 0.0001 inches, and if the pressure of 


| ane mer 


the gas stream is approximately 100 microns, the ratio of the mean free path 
F to the wire diameter will ^e approximately 180. Therefore the flow around the 

wire is definitely the free molecule flow (Ref. 2). We have thus a simple phy- 

sical situation, which is an improvement over the rather uncertain circumstances 

of mixed dynamic and visoous effects for the measurement of velocity by Pitot 
tube, It thus seems worthwhile to explore this possibility by a trial calcula- À 
tion of the performance of such a hot-wire anemometer, 
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& macroscopic velocity U and а Maxwellian molecular velocity distribution, 
the translational energy of molecules incident upon the unit ares 16, 


(аз) 


where c^ = 2 RT, T temperature of tho gas stream and erf is the error function. 
Now let г be tho radius of the hot-wire. Then the total energy B incident up- 
on а unit length of the wire ів the sum of translational energy and internal 

energy. If Q, is the specific heat at constant volume, this total energy per 
unit length of wire is 


(a9) 


The integrals in equation (19) сап be expressed in terms of tabulated functions, 
(see Appendix) Thus 


(20) 


(21) 


(22) | 


The I, and I, are the modified Bessel functions of the first kind of ordere sero 
and one respectively. The functions P, and Р. are tabulated in Table 1. 
If T, is the wall temperature, and x the accomodation coefficient, the 


difference „СОЈ. МЕ and the energy 


Ἐς carried by the molecules re-emitted from the surface is given by 
қ-к- (N - E) 


whore By is the energy that mould be carried away by the molecules if the 
re-emission мего at the temperature Т, of the wire. Therefore 


қ-қе cr dU. Gre ( 20 f (2) 


+ aveo a- nb а» 


This difference of energy is then the net energy input to the wire per unit 
length of tho wire by the air stream, 

If i is the electric current heating the wire and is the resistance 
of the wire per unit length at the wire temperature, the heat input per unit 
length of wire by the heating current is 1^ . Hent is lost from the wire 
by radiation. If is Stefan-Boltsmann constant and “is the emissivity of 
the wire surface, the radiation hent loss per unit length is š 
Thereforo if the wire has reached а steady condition, the heat balanoe re- 
quires 


(24) 
This equation can be put into somewhat simpler form by using tne relation that 


= c- . )په‎ -1) (25) 


OA nnti 


% “nê the corresponding temperature coefficient of the re- 
sistance is . Then the resistance олп be expressed as 


= 01» (%-%,) (26) 
(27) 


а +1 (28) 


o 
Mow introduce „ ав the reference density and 1, as the reference heating 
current, then equation (24) oan be written аз 


(29) 


‘The partiouler values of the reference temperature То, the reference den- 
sity ϱ and the reference current i, are not yet fixed. We fix these quanti- 
ties now by requiring that 


“41 (30) 


(a) 


а ОРЛА 


(33) 


This is then the performance equation of the hot-wire in free molecular flow. 
Now let us investigate in greater detail the oase of a bright platinum 
wire. To satisfy equation (30), 


To = 492? n. 
The value for and сап be taken to be 0.06 and 0.90 respectively. Then 
equation (31) gives the corresponding pressure p, for ο and T, ав 


Let the radius г of the wire be 0.0001 inch. Then equation (32) gives the 


reference heating current i, as 


14 ------ = 027 silliampere. 

where the resistivity of the platinum is taken as 10.96 х 1076 ohas-oas. There- 
fore the order of magnitude of the different quantities 18 entirely satisfactory. 

If the wire is used with a constant heating current, then equation (33) can 
be used to calculate the relation between the resistance ratio and the velocity 
төне ( B.) st constant air өзгені density and tenperaters, This 19 done for 

RS a and the result is given in Fig. 6. It is seen 
that the sensitivity of the instrument is good. Of course, the behavior of the 
hot-wire anemometer will be actually determined by calibration for any experiment. 


ж ‘The author is indebted to Mr. L. Mack for the numerical computations. 
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for the model test to be similar to the prototype, the model must be made of 


same surface material as the prototype, and the fluid must be the same, and 


furthermore tho following parameters must be the same: 
(1) Reynolds number Re 
(2) Mach number w^ 

(3) free stream temperature T° 
The radiation heat loss from the surface is equal to 


per unit area, 


However, if the model is surrounded by the walls of the test chamber, there is 
also an heat input due to radiation from walls of the test section to the model, 


Let us call this quantity . Then the net heat loss per unit area of the sur- 
fnoe of the model is + This quantity can be rendered non-dimensional 
by dividing it by + Call this new parameter , then i 


(34) | 


For the prototype, the heat from the walls of the test chamber is absent but 
there may be solar radiation and the radiation from the earth and surrounding 
atmosphere. (Ref. 0). Denoting this amount by , then the parameter Гог 
the prototype is 


(35) 
In order for the flow to be same also with respect to the radiation heat transfer, E 


(36) 


— 


Because of the previous conditions on the Reynolds number and free stream ten- 
perature, (36) is the same as 


te o 


coefficient аз shown by equation (29), it will be necessary to find mater- 
ials which can hold this coefficient constant for a considerable period of time 
во that no frequent calibration is required, However the present analysis seems 
to indicate the feasibility of such an instrument for measurements in rarefied 
gases and further research is definitely desirable. 


4. Parameters of Flow 
The two parameters that are directly connected with the flow field 


are the Reynold's number Re, defined as 


nat 


where 18 the kinematic viscosity, L the typical linear dimension of the body; I 
and the Mach number M^ of the free stream, “This is true even for slip flows and 
free molecule flows due to the fact that the ratio of mean free path to the typi- 
oal dimension can be also expressed in terms of tho Reynold's number and the Mach 
number. 
However, аз the pressure or density is reduced, the solid boundary of the 
flow enters actively into the flow oonditions by requiring not only that the 
microscopic stream velocity be tangential to the surface but that the interaction 
of the molecules and the wall be considered and that the radiation of energy to 
and from the wall be taken into account. Тһе interaction of the molecules with 
the wall is so far expressed through the fraotion of molecules that are dif- 
fusely re-emitted from the wall, and the accommodation coefficient . Itis Т 
known that both and аге functions of the temperature of the wall and we have 
reasons to believe that they are also functions of the molecular velocity distri- 
Therefore the interaction of the molecules with the wall iə the same 


bution. 


only if the wall temperature, the gas temperature and the Mach number of the gas 
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The Functions P) and E, ( of. equation 21) 


We қ 


СОРҮ 


Braluation of the Functions P) and Р, 


Yor the funotion F), 


whore I, is the modified Bessel function of first kind and under sero. Тһе 


last step із made possible by the substitution LI 
For the function Рр, 


By the definition of the error function, we have 


‘This form oan be simplified by partial integration. Thus 


where Ly EUR Q | . 


linear dimension of the prototype. This means that the wall tenpera ture of 


ths test chamber must be so controlled that satisfies (37). 


This set of ratuer strict similarly rules for model testing in super- 
aerodynamic flows is certainly difficult to satisfy, In what way the rules 
can be relaxed is the problem of future resvarch. 
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Therefore 


where Ij is the modified Bessel function of first kind and order one. 
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The problems in the experimentation of rarefied gas are discussed. 
First the extremely large viscous effecta іп а wind-tunnel nossle is 


shown. Then the difficulties of flow measurement are surveyed, pointing 


out particularly the unconventional behavior of the Pitot tube in 
rarefied gas. The performance of а hot-wire anemometer is then studied 
in some detail to show its feasibility. Finally the rules for achieving 
complete flow similarity of rarefied gas flow are formulated. 
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Application of Tschapligins Transformation 


to Two Dimensional Subsonic Flow 
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Section 5 


The Buckling of Thin Cylindrical 


Shells under Axial Compression 


THE BUCKLING ΝΌΒΙΟΛΙ, SHE] UNDER 


Th. von Karman and Hsue-shen Tsien 
еме Adi 1% 

In two previous pepers (Ref.1 and Ref.2), the authors 
have discussed in detail the inadequacy of the classical 
theory of thin shells in explaining the buckling phenomenon 
of eee cylindrical simi end the spherical shells. 

ТЕ was shown that not only the calculated buckling load ia 

4 to 5 times higher than that experimentally observed but 

the buckling wave pattern found is also different from that 

predicted. thet the 

different explain&tions for this discrepancy advanced by 

ын. Donne! P u. Flügge ee untenable when certein 

conclusions drawn from these explainations are compared with 

the experimental facts. — The authors are then led by 

bash ta theoretical investigation on spherical shells (Ref.1) 
а model experiment (Ref.2) 6n “thin eólumns 

with non-linear elestic support/to the belief that the 

Duckling phenomenon of curved shells can only be explained mv 


moans of tam non-linesr zarge deflection theory. 
tae laite κ. 


to keep the shell in equilibri: 
buckle. This chsrscteristic of „ load shows, first 
of all, that there is a release of elastic μα ον the 


buckling has started, and thus explains the observed rapidity 


2 
Furthermore, this characteristic 
also brings in the possibility of a decrease in, buckling 

load when there e slight imperfections in the test specimen 


of the buckling process. 


te nen there ёге vibrations during the testing process. 


In this paper, the authors will show by means of 


an approximate calculation the dropping load characteristic 
ние of a uniform thin cylindrical shell under axial 
compression. Consequently, they hope that they have thus 


offered ап aceptable explaination of the observed facts. 


Stresses in the Median Surface and the 
Expression for the Total Energy of the 
System 
(39) 
τε κ, and. тұраты ἀῑκρταζφώνος DE a polnb-ón Mu 
median surface of the shell in the axial, ху direction, 


the circumferential, yg-direction, and the radial direction, 
then the unit strains in the x snd y directions, fy, &y 
and the unit shear ty at a point in the пөйійп surface 
сеп be expressed іп the following form, including terms 


up to ke δα е T 
%- FH 5 Ht er be) 


where R is the radius of the undeformed median surface of the shal. 


The strüses an tre ins in the median surface of the shell 
аге, however, related to each other by NL ons 


= uet s v < aina, pats, 
eats pan асы Eq.s(1) into Equs.(2) the following 


μα bse between the Stress n the median surface 


"t tne `S Кене of the median surface are obtained: 


=] Ë 7695 "t 9 M-E | R 
2 . ο жі | 
τε is 7 / ἡ Helly accepted 4/6 the equilibrium conditions 
of the stresses іп the median surface can be RS Sg 
expressed by the following equations in case of thin shells: 


mt EN fw 


(4 
7 А я; ға 
This pair of equations KA be амімеміні-еефіу satisfied by 


introducing the well-known Airy stress function, F, defined 
by the relations 


a og pd = 


Eliminating the variables u and v between Eos. (5) and (5) 

lation between the Airy stress function,F and the radial алайы 
a rela’ pop ry E ang, 

ас! з obtained; 


ΤΝ 42,53 _ Z 22 
25. 5%. dre e P] FEF: {δω 


4 


This equation 18 very similar to the large deflection 
equation for flat plate derived by the senior author. ( 49 
L.H.Donnell (Ref; 3) first obtained this equation ip its 
EPS 
present form. With a given form of the radiki қ 
Eq¥.(6) gives the induced stresses in the medjer surface of 


uon 2421 же e egt e aon له‎ the 
vaca. To epson 222 elastic energy, it is 


төс! 


гу first to find the expressions for the change of 


or curvatures % twist of the median surfaci H 
П In this paper, the following simplified expressions зей: 
24 Ya y 
AC = 22 - X 24. 
b t ap, % ga у 27 (2) 


τη Eq. (8), certain additional tems in Yyand X., involving 
v are neglected. | It was shown by L.H.Donnell (Rer. B) 

that those tams b... differ from the terms retained 
in Eq.(8) by a factor Jf} where п is the number cf waves in 
the circumferential direction., For thin cylindrical shells, 
the value of n is r around 10; therd the neflection 
їз jusitfied, With these expressions for the change of 


the curvatures t d the median surface, the bending 


fe ее кий те сі 一 


Gan be obtained as 


Бей 9) Бага v 


L^ 
, 
The lowering of the potential, of the applied force 

on the end of the cylindrical shell can be calculated as 


the product of the change іп length of the shell and the 


applied force. 22 a ғеу 


k en - 


κα: 7/8 E s 


The equilibrium condition of the shell із obtained 
by equating the first variation of the sum of &néPgies, Wy Wo 
and Ws to zero. Іп the following calculation, the Rayleigh- 
Ritz method or the energy method 8813 be used by assuming а 
plausible form for the function ж. 


Calculation of the Total Energy 


To obtain a plausible form for w , one has to resort to 


the Tesults Swemlexperimante$ It 1з observed that, for large 


values of the wave amplitude, the waves are so-called diamond 


shaped. THis particular wave shape can be approximately 


πω. ба бт - 14) Liu) 


σα 


where the squares аге introduced to account for the fact that 


expresséd by 


the shell has a definite preference to buckle inward. Eq. (п) 


6 
сап be rewritten аз 


бетер nt e o 


On the other hand, the classical theory which is correct for 
infinitesimal velues of the wave amplitude, requires the 
waves to be of tne form 
M qu, DE 
lo 4 
e e ασ) 
In order to satisfy this requirement, the wave form assumed 
in the following calculations 1s of the 


-- che- [n Же f. fe fe 
+ FA [ta + en E] 
where fo, fi, ата fa are unknowns to be determined by 
conditions of minimum totel energjes of the system. 


fj is introduced imfrder to allow the shell to expand 


o мый Аё же tak а- la 56 44 acent 
radially. Ж 2 Aid — p. cerros 
7014) into Eq.(6), the differentia 


equation for the Airy stress function F is obtained: 
4 B cu Е 
лак gl e, Hn n 
+E ane e Ж + C cat ta Ж +F er , و‎ 


: 


ναι uts g É — 


/ tl 


The solution of this differential equation can be easily 
obtained as 


4 М 
эё * gz 
ge- FH P m uM. m ТАРДА 
+Z. P жаз 7:72 n TR ngu, pre 
Using Eq. (5), the stress components in Б median ML 


— аз о. de ὧν 
%- ep Боч uen фаре 
+ Бремен, ep %, й 
ар ма 
σ- πα ЖР to tn + η” ro 
2 эм ey? 
ie екеу dis ento < 7 
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In all experimental results, the deta із usually 


expressed intterm of the average compression stress in the 


axial direction. It сап be easily seen from Eq.(18) that 


"T (19) 


Вғ-есіліпд Using Eq.(3), t dor the following 
expressions for 25 ona BF can be obtained: 


ΙΝ 
%- Ef v -F K. 


By substituting Eds. (14) and (18) into Fa. (20), it 15 found 
that 


ж $ > 2 tine Term. 
y - dpt bli neu] mme 


(21) 
у * Жет S -fn "ж pt pelr H) 4 rotine Term, 


Since -en (ее circumferential direction, v 


must be а periodic function; therefore the constent te m in 


mat be Ы Sem Or 
p 
E v£- /I V ALIE 


This ta the сопан tone’ aeterninele foc. » 


c 
Using Eqs. 416) and (гә), the externsional energy of 
the shell is obtained as 


dix bed, "н 
ο.) e Ax terre, 

2 D “ ss # * $ $ 
EE an ffir pe 


“ғ 


Using Eqs. (9) and (18),,the bending energy of the shell 
can be calculated as 


Ж pa У (14 
πα we I pP 
= Op th + (rh | (24) 
The lowering in potential of the afl 1c force can be 


obtained by means of Eqs.(10), (18),£421}2ada (22): rh 
result is 


Za ae [rime а pent P 
pnt үрү юе nei 


Relation between the Compression Stress 
and the litude of Waves 


To find the relation between the compression stress 
and the amplitude ФГ the waves, the conditions which will 
make the sum of the energies W, ,Wo , and Ws в minimum 
have to be obtained. It was found that the calculations 


can be simplified to в certain extend by first жә = а A. 


көн (28) 


This condition тір determine 5 relation between fo 
ana f, end fg y 


It vide Ë т 


Subdituting Eq.(27) into the expressions Гог W, ,W р and Ys 
as given by Eqs. (25), (24) and (25) and using Eq. (16), 

the total energy Ж the system 15 expressed finally in 

the following form: 


NM gr gyal 1 4) 
ως * is * 

2 tl ( 8 * 4 Z t mH 

(ER — ZJ ll bert eF GI 


& 
ЕСЕ Е 
а dup + др d E 


11 
The equilibrium conditions are then obtained by differentiating 
the expression for té&al energy, Eq.(28), with repppct to 
Ті end fj , and then set these expressions equal to zero. 
The results can be written in a simpler form by introducing 
the following parameters 


- 4/4 

j= ”# (29) 

ЕЕ 
人 


Tis ше erii oen edn are 


3E μή.) - ife А йу” СЫ 


12 
Eliminsting τή ffom the above equations, the following 
equation for ic is obtained 


Apt Age Apt 452 


where the coefficients are 


ssl 


Αρ 


(52) 


Thus with a given value of / and Z , the coefficients 
for various values of the wave amplitude ғ сап be ftrst 
calculated by using the Eq.(352). Then Eq.(31) сап be solved 


for different values of [4 corresponding to this particular 
value of Жала / at vartous wave amplitude $ . when 
the value of Ф is known, Eq.(30) can be used to calculate 
the corresponding value of the compression stress Е. . 
It 1s found, however, that the following expression for the 
compression stress which is obtained from Eq. (30) by 


15 
‘eliminating the third powers òf 3 „18 more suitable for 


numericel computftion: 


These calculation was carried out for two values of 
the parameter 1 — —on—ono 
бөте the ratio of the wave lengthscin circumferential 
direction and in axial direction. These values of / are 
1 and 0.5. The former value was chosen, because the 
experiments indicates that at large wave amplitudes, the diamond 


waves hat qual sides. The latter value of JL was chosen, 
because the experiments made by N. Nojima and S.Kenemitsu 
(Ref.2) shown the occurance of tnis type of waves. The 
results of these computations àre shown in Fig.2 and Fig.3 


where the compression stress is plotted against the wave amplitudes, 


14 
The parameter in the figures is 

/ we values written 
in the parenthesis is the actual number of алыс; circum- ` 
ferential direction fer - = 1000. Тһе dropping 
load руе is VAI evident . 


474. 


the actual losd taken by the specimen under a testing machine 
can not be predicted from these figures. In a testing 
machine, the only fector under the control of the operator 

1s the distance between the end plates and this is the feeeor | 
whishidetwsmine geometrical restraining the specimen 

must confgÉm. Therefore, it is clear that the compression 
stress on the shell 1s not determined by the wave amplitude, 
but by the end shortening. The unit end воен 1.е., 


from Eq.(21). It is found that 


f = Bs НИР Ф) 


Ж. στην 
Ж fe- РУС XC 
arial thene oi ah Ri А. 


15 
Thus the unit end shortening canbe easily calculated by 


4 and 5, the compression stress is plotted кеа the 
unit end shortenings, for . / апа 
It is immediately clear from this two figures that if 
the buckling process follows the curves iram once the 
shell starts to buckle the erro or tege 5 
the end КЕТЕ to ae Art. Thais were 
Ps ^ev rode e nin: peter 
— £ ee te „4 A 


‘Tharaforer the shell must £¿rst jump to the point P (Fig4 or 
5) which has the same end shortening as the starting point of 
buckling process, but has a much lower compression stress . 
This jump definitely involves в release bf elastic energy 
and thus explains the rapidity of buckling process observed 


——— ae fX accompanied vibrationy zb же afte 2 
However, it is still difficult to ss 9 shell will 


jump Mo в state with narrow waves, as shown by Fig.5, аг 
Οτο а state with almost square waves, as shown by Fig.4^ 


Superficially, one might conclude that narrow waves е 
probablė because it gives a πο e618 stress, аз 

can be seen by comparing Fig.4 with Fig.5. But at a 

given end shortening, [the true criteri 

fundamental consent of mechanics/for the most probable 
equilibrium position 18 that the elastic energy stored іп 
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shell should be the lowest. An арргохітьбе calculation 


of the elastic energy of the shells shows that for values of 
/< 2.0121 » the occurance of narrow waves is improbable, 
because the elastic energy is higher than that [Z square waves 
at same value of end shortening. | However, for fÉ near 
0.6, the elastic energy stored 1n the shell for narrow waves 


1s comparable to that for the square waves at same values of 

the end shortening. This indicates the ssibility of fh ana 
мо intial ж ا‎ Кине, 1 

narrow waves 

— | | This is in agreemnt 

with the experimental findings of N.Nojima and S.Kanemitsu 


ав reporte? in а previous ра) by the present authors (Ref.2). 


In any οἱ 
positions of a buckled cylindrical ммм which involve 


it 1s defámte that there are equilibrium 


much lower average compression Жейт се that sroga 
— owe 

the beginning of buckling. Teke for instance, the sqaure 

ivo a lowest compression stress 18 given by 


2 р-н 5 


This value corresponds closely to most 1 experimental 
results obtained by E,E.Lundquis т.н. Donnell, (%) 
In the two previous papers (Ref.l und Ref.2) the authors 


have pointed out that due to imperfections present in 
test specimen and vibrations in the surroundings, siats 
preebles+ the test buckling load will probably nearer to 


17 
the minimum load necessary to keep the shell іп buckled 
shape rather than tne Wr buckling load given by tne 


pes gee theor 


The Effect of the Elastic Characte istic of the 
Testing Machine on the Buckling Behavior 


In the previous paragraphs, it 1s stated that the 
compression on the specimen 1s determined by the distance 
between the end plates which is in turn determined by 

the loading mechanism. However, in actual cases, thore is 
always а certain amount of elasticity in the testing machine 
and thus the distance between the end plates is not only 
determined by the loading mechanism but also by the compression 
force weh on the specimen. With the loading mechanism 


E. at в certain position, a compression forge on on tne specimen 


will force the eng plates spart and thus төле е Ond ‘shortening 
hed. < Pa — 
Jer the Ape τα that the testing ma tj 
а linear elastic behavior, the &ompression load is related 
to the end shortening by parallel straight Inez Hör Constant 
values of, s&g,-number of turns of the loading crank of the 
testing machine. If the léading crank is held at a fixed 
position, then the compression load of the specimen must 


follow the straight line corresponding to this crank position, 


18 
Іп Fig.6, the load-end shortening curves of different 


character are superimposed on a system of straight lines representing 
fhe characteristics of 


ting machine. It is evident that 
after the maximum, or initial buckling load is reached, the 
shell will jump to ax cutlibriun state involving much lower 


compression load. e new equilibrium state 1s determined 
by two factors: ΜΕΥ of the load-end shortening 
curve, Тһе curve А la pre load than 

the curve B (Fig.6). elastic character of 


the testing machine. А more elastic machine will give 


а set of characteristic straight lines having smaller slop 


Therefore, in ca: 


of Curve A, a more elastic machine will 
result in а higher load, while in case of curve В, a more 


elastic machine will result in в lower load. Therefore, 
Jt 

the buckling behavior of a cylindrical shell 1376117 deterhined 

by the elastic characteristics of the shell itself, but also «а-а. 


by the elastic characteristics of the testing machine. | 
| 
| 
| 


< өзенде — 
This fact is first discussed by the senior author imis 


dut eec (Ref. b. 


Tftheve“is-e-slight imperfection tn thetesting~speeinen " 


then 
Concluding Remarks 


у Іп the previous paragraphs, the authors have shown 
that there are equilibrium positions in the buckled shape 
involving much lower load then the buckling load predicted by 
classical thoeyy, and thus if the specimen is slightly 

imperfect, it 15 reasonable to expect much lower buckle loado 


quA 


ma 


ord ved in the laboratory. They have 


also pointed out that the elastic behavior of the testing 
machine has a profound influence on the buckling phenomenon 
and this might be another cause of the large scattering of 
E obtained by different experimenters. 
However, due to complexity of the problem, the results 

given in this paper can be only considered as a rough 
μμ eei οὔ Ge айына wade Se μια 
rather than quantitative. To put the new theory on a 
solid footing, a more accurate solution of the differential 
ecuations of equilibrium is nei агу. Particular attention 
must be given to the càlculation of the elastic energy stored 


in the shell, because it is found ао — 
. ee σσ 


. Г the compression stress on the specimen nor the wave amplitude of 


the buckle? cylindrical shell, but rather the elastic energy 


stored: ыл < 


< 
Rerheps, а тиме inquisitive trreettgetor е pleased 
by а rigorous proof of the validity of all the large deflection 


equations. equations weed are 

g intuitive arguments rather by systematic rı 

Рог instance, it is not certain whether the curvature of the 
shell has to be calculated more accurately by taking into account 
the second order terms, or the extensions of median surface 
should be more accurately determined. | It is the belief of 

the pmeeert authors that an investigation of these problems 


20 
by starting from the general non-linear elasticity theory 
бозаега® (Fer.9) and others e very Haney Даа 
The sentor author has already expressed this get 
opinion in his 1949 01005 Lecture A. 1 £ du L.E t 


bak 


TEE BUCKLING OF THIN CYLINDRICAL SHELLS 
UNDER AXIAL COMPRESSION 


Theodore von Kármán and Hsue-shen Tsien 
California Institute of Technology 


In two previous papers (Ref, 1 and Ref. 2), the authors have discussed 
39 detail the inadequsoy of the classical theory of thin shells in explain- 
ing the buokling phenomenon of cylindrical and spherical shells. It me 
shown that not only the calculated buckling load is 5 to 5 times higher 
than that found by experiments, but the observed wave pattern of the 
buckled shell is also different from that predicted, Furthermore, it was 
pointed out that the different explanations for this discrepanoy advanced 
Әу L. H. Donnell (Ref. 5) and W. Flügge (Ref. 4) are untenable when certain 
conclusions dram from these explanations are compared with the experimental 


facte, He rattrore -are -then lot By a theoretical investigation on spherical 
уриа рут 


shells (Ref. 1)/to Um pedet 3 of curved shells 


oan only be нвн зу means /of non-linear large deflection theory. This 
point of view tr substantiated by 4 model experimentson “rin vorm 
with non-linear elastic support (Ref. 2). ais ады ы ا‎ aE 
Anvestigatione—thet the load necessary to keep the shell in equilibrium 

drops very rapidly with increase in wave amplitude once the shell started 

to buckle, This characteristic shows, first of all, that there is & release 

of the elastic energy stored in the shell once the buckling has started, 

and thus explains the observed rapidity of the buckling process. Furthermore, 

4% also brings in the possibility of а decrease in the buckling load when 

there are slight imperfections in the test specimen and when there are 

vibrations during the testing process. 


In this paper, the authors will show by means of an approximte 


-ᾱ- 
calculation that in case of а thin uniform cylindrical shell under axial 
compression, the load sustained by the shell drops with increasing deflection. 
Consequently, they hope that they have thus offered an acceptable explanation 
of the observed facts. 


Stresses in the Median Surface and the 
Expression for the Total Energy of the System 
Let X and 4 be measured in the axial and the circumferential direction 
їп the median surface of the undeformed cylindrical shell and /, V , and 


W be the components of displacement of а point on the median surface of 
the shell in the Х direction, the f ~direction, and the radial direction. 
Then the unit strains in the X and нения, & & % and the unit 

shear. Ú, at в point in the median surface сап be expressed in the following 


orms, io luc iE ga + 
d я 5 1p € а | 


e HEE 262 A 

where © is the Jn of the undeformed median surface of the shell. The қ 
stresses and the strains in the median surface of the shell are, however, 

"4 related to each other by the following equations: 

Q= Er (a 

Жады Ашы e 


Ty = zm бу 


b 
where E is Young's modulus of elasticity and ¥ is Poisson's ratio, There- L 
fore by substituting Eq. (1) into Eq. (2), the following connections between 4 
‘the components of stress in the median surface and the components of die- 
placement of the median surface are obtained: 


6-2. Hl: ty #/ 
ер" rope Xl 
κ σε, ЖЕ HH] 


It is generally es that the EL conditions of the 
stresses in the median surface of & thin shell oan be approximately 
expressed by the following equations: 
20 А 254.) 
3% + 3 «ο 
This pair of equations can bé satisfied by introducing the well-known 
Airy's stress function, Fle f) ажалға by the relations 
τ N za s= 92; 
E, yy Fo 
Eliminating the variables Д and 7 in Eqs. (5) and (5) a relation between 


Airy's stress function 2/Х, p= the radial component of the displace- 
ment, W is obtained: 


ο τ “ERD vy ΚΕΙ в) 


This equation ο ‘the condition of compatibility between stress and 


(4) 


strain, When f— «ο, it reduces to the corresponding equation for а flat 
plate derived by the senior author (Ref. 5). L. H. Donnell (Ref. 3) first 
obtained this equation in ite present form. With а given form of the 
radial component of the displacement, /, Eq. (8) gives the induced stresses 
in the median surface of the shell. 

For one complete wave panel, the extensional elastic energy LA corre- 
sponding to these stresses can be written as 


44 
жі ολ 3/% 00 


-t 


where 2 «nd / are the half wave lengths in the axial and the cirounfer- 
ential directions respectively. 


To calculate the bending elastio energy, it is necessary to find the 


expressions for the change of curvatures and the unit twist of the median 


surface, 


In this paper, the following simplified pum will be used: 


2 ge ‚ №- ЖУ (в) 


In Eq. (8), certain additional м LA Жу ыр are 
neglected. It was shown by L. Е. Donnell 24 


6) that those neglected 
terms differ from the terms retained in Eq. (8) by a factor = „ where 
JL is the number of waves in the circumferential direction. For thin 

cylindrical shells, the value of 7] is around 10; therefore the neglection 


is justified, With these expressions for the change of curvatures and 
the mit twist of the median surface, the bending energy ДУ for one 
complete wave Т” can be written as 


Wait, 7] ait, ή / Жж en) Hj VE Bf (ө) 


V 
of the cylindrical shell can be calculated as the product of the change in 
length of the shell and the applied force, Therefore the following 
expression is obtained for one complete wave panel: 


W= -14 αι 4/44 loj 
Tho equilfbrium condition of the shell een be obtained either by 
equating the first variation of the sum of the energies LAE жа Hy to 
zero, or by actually analyzing the moments and the stresses in the median 
surface of the shell. Using the approximations stated previously, Donnell 
(Ref. 8) derived the equilibrium equation ав 
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A i 2 „ 
Zen lite м Aach. 您 ЖИН 24% 
where fj is the external radial pressure оп the surface of the shell. In 
the case omcerned, $ =0 » then using Eq. (5), the second equation 


connecting Airy's stress funotion FÉ f Joos the radial component of dis- 
placement, A/ is obtained as 


Жө eed d ЖЖ %% ЛЕ 


‘When еб; Eq. (12) reduces to the corresponding equation for а flat 
plate. 

There are two different ways to solve the problem of the buckling of 
A thin uniform cylindrical shell under axial compression. The more exact 
method is to solve Eqs, (6) and (12) simultaneously, using appropriate 
boundary conditions. The approximate method is to first assume a plausible 
function for J , with mdetermined parameters and then use Eq. (6) to 
determine the stresses in the median surface of tho shell. The energies 
И. W, saa № сап then be caloulated by means of Eqs. (7), (9), and 
(10). The undetermined parameters сап be ascertained by the condition that 


the sum of energies M; , N. Ж) must be a minimum, This approximate 
method will be used in the following oalouletions. 


Caloulation of the Total Energy 
To obtain а plausible form for W , one has to resort to the experi- 


mental results. It is observed that, for large values of the wave amplitude, 
the waves show а so-called diamond shaped pattern. This particular wave 
shape can be approximately expressed by 

ML ume (men 

κ JE “κ 
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where the squares are introduced to account for the fact that the shell 


has а definite preference to buckle inward. Eq. (15) can be re-written as 


. lle keln. de αν 


Оп the other hand, the classical theory which is correct for infinitesimal 
values of the wave amplitude, requires the wave to be of the form 


* ¿e TE ou 2E (ав) 


In order to satisfy this requirement, the wave form assumed in the following 
calculation is 


* 2 Жетен poet Шет T) 


were 7,24, μος 

minimus total energy of the system, 4 is introduced in order to allow 
the shell to expend radially. Тһе amplitude of the wave pattern defined 
as the maximum difference in the radial deflection, J/ is evidently given 
зу 7. The wave lengths in the axial end the circunferential direction 


are A. and zz respectively. Hence the number of waves along the 
circumference of the shell is equal to y . It is evident that no end 


effect can be accounted for by this form of wave pattern, and therefore 
‘the following calculation really corresponds to the oase of a very long 
cylindrical shell, This simplification is justified by the experimental 
findings of N. Nojima and 8. Kanemitsu as reported in a previous paper 
(Ref, 2). It was found that there is no appreciable length effect when 
the length of the cylindrical is greater than 1,5 times the radius of the 
shell. Furthermore, it is seen that by setting f . Eq. (16) is 
reduced to Eq. (14); эпе зу setting Z, Zp ea 4Е- „ Eq. 


(16) is redueed to Eq. (15). With other values of these parameters, wave 
patterns intermediate between these two limits can be obtained. 


Substituting Eq. (16) into Eq. (6), the differential equation for 
Airy's stress function ғу) is obtained: 


(GERNE - gell hene В cadet Соо сед 
+0 en len By Gen Ben +H cn ge ., = 
where A A. » the aspect ratio of the es. If > / the πάτον are 
longer in the circumferential direction; it <1 » the waves are longer 
in the axial direction. The coefficients in Eq. (17) are given by the 
following relations: 
as phn- HH) 
B= gn 
oh 
s" An 444) 
ο- {4144 +h) 
se EE 


The solution of this differential equation can be easily obtained as 

ΕΕ Бора tat i + фе entea% 

лет gio nnt жі 19) 
speedy 

Using Eq. (5), the stress components in the median surface can be written 

Келелі meal int 

MM MAUS 


(20) 


5- a d “желеу айтад 
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the average compression stress Û in the axial direction. It can be easily 
seen from Fq. (20) that 


Bre. (a) 


Using Ра. (3), the following expressions for Жж and LE oan be 


obtained: Й 2 2 
-Z -A, 
2 (22) 
Lp ve ЖО" 
By subst} Eqs. (16) and (20) into Bq. (22), it is found that 
4 
B [era n HAL ΚΟ} mt eria 
т 5 (23) 
ж fef £v: gat E Terns of periodie 
Since / ів measured along the circumference of the shell, И must be a 
periodic function of / 1 therefore, the constant term in фе must be 
equal to zero. Or Z 
ЁНИ но ш 
This is the condition for the determination of Æ . 


Using Eqs. (7), (20) and (24), the extensional energy Же the shell 
is obtained as 


ED Ξ ЕТТЕ И (NE at 2 rh bf 
ae BEA ASMA f+ ΤῊΝ 
Ty e Ep.) 


Using Eqe. (9) and (16), the bending energy //, of the shell oen be 
oaloulated as 


жа” ο нее ЖОЛУ 


of the applied force сап be obtained by 
means of Eqs. (10), (20), (23) and (24). The result is 


κα gp rtp 
£u th) ] 


Relation between the Compression Stress and 
the Amplitude of Waves 
То find the relation between the compression stress and the amplitude 
of the waves, the conditions which will make the sum of the energies ,. 
И; + ena И а minimum have to be obtained, It was found that the 
calculations oan be simplified to а certain extent by first using the 
condition that the sum of energies must be minimum with respect to . 


Е (тен ж) -о s 


This condition determines a relation between A and £ and 4 , which 


can be written as: 


4% - ales tth #47 (29) 


Using this relation and Eq. (24), it is easily seen that 4/72 „ In other 
words, the shell will expand radially to such an extent that the average 

of the oircumferential stress 9 is equal to zero. Substituting Eq. (29) 
into the expressions for Ж) , Ж and Wy as given by Ede. (25), (26) and 
(27) ап4 using Eq. (18), the total energy of the system is expressed finally 
in the following forn: 


i UEP- Fir (Ph AA “4614755. 


+ dates td ΚΩ 750 a η xd #5) ud 
dies фи К ge, АСА Z dept] 
- п | Le πως 
dL ird тия fs F. (во) 
АЛЫПТАРЫ ЛАРА РАУ) 


The equilibrium conditions are then obtained by differentiating the 
expression for total energy, Eq. (50), with respect to f end 4 , and 
then set these expressions equal to zero. The results oan be written in 
а simpler form by introducing the following parameters: 


rior т en 


" 


where Š is the wave amplitude of the buckled shape of the cylindrical 
shell. Then the equilibrium conditions are 


Жез ту er e, 
[tr , Hi, tish lE l E, ts tiol 
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Blininating OF fron the above equations, the following equation for f 
is obtained 
3 
ASHAR -0 


(55) 


where the coefficients are 


ws Oy fes “5, 4, 
V фе, ЙА ЖАЛА 
0% l be, "б e ЦЕ ΑΡΗ 
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Therefore, when f =0 , i.e., when the wave amplitude approaches rero 
ва. (82) gives 4-4-0 1 and 


4- ης 2 
κα. 


Substituting into Eq. (31) it ie seen that -g. vr fats. 
Putting this relation between / snd {into Bq. (14), the wave pattern is 
reduced to that represented by Eq. (15), i.e., the classical wave pattern 
for infinitesimal wave amplitude. 

With a given value of Ж and 7 , the coefficients for various values 
of the wave amplitude Z oan be first calculated by using the Eq. (34). 
Then Eq. (35) can be solved for f? corresponding to this particular set of 
values of f and 7 at various wave amplitude f . When the value of 
Ав known, Eq. (52) can be used to calculate the corresponding value of the 
compression stress g + It is found, however, that the following 
expression for the compression stress Ж which is obtained from Eq. (30) 
by eliminating the third powers of e is more suitable for numerical 
computations: 


ЭА 


H- [tir to A 
vello ση Же) E 
Coal tay tite doge] ONE tag 
enti н Alte: 
ҮЛЕМ ЛІ 


As- 
Therefore, when F20 , i.e., when the wave amplitude becomes very small, 
Eq. (35) reduces to 


А 
— 
pat 

ЕР „=, 7 ff, ge, jf (зв) 

The minimum value of T average τ stress is given by 
№ | Z „ ТҮГІ; Be (37) 

which із the well-known rine from the classical theory of infinitesimal 
deflections. This minimm value is obtained when 


2 
ر‎ = M 30v] вв) 

It is interesting to notice that for infinitesimal wave amplitude, the 
minimus value of average compression stress is not determined by separate 
parameters 7 нд „ but by a combined parameter shom in Eq. (38). 

These calculations were carried out for two values of the parameter 
JM. » the ratio of the wave lengths in cirounferential direction and in 
axial direction, These values of М are / ana 25 . The value / 
was chosen because the experiments indicate that at large values of wave 
amplitude, the diamond waves have almost equal sides. Тһе value 25 
ms chosen to investigate the possibility of occurrence of narrow waves. 


The results of these computations are shown in Fig. 2 and Fig. 5, where 


the compression stress 25 is plotted against the wave amplitudes . 


The parameter in the figures is / . The values written in the parenthesis 
after ) is the actual number of waves % in cirounferential direction for 
R/E = 1000. Fora given value of 7042, i.o., а fixed size of tho wave, 
the load sustained by the shell, g first decreases as the wave amplitude, 
E , is increased, After а minimum is reached, the load will rise with 
increase in wave amplitude, When the waves are larger, the initial buckling 


Er 

load, i.e., the value of 2f at Z= , is higher. However, the minimum 
load reached tends to a lower value, except for /<4//Уала 4/0. For 
А05 , the lowest value of the minimum load is not yet reached at He c- 


The Relation between the Compression Stress 
and the Shortening in the Axial Direction 


Although the load characteristic of the cylindrical shell shown in 
the Figs. 2 and 3 give the possible equilibrium relations between load and 
amplitude of the deflection wave, the actual behaviour of а specimen in a 


testing machine cannot be directly seen from these figures. In & testing 
machine, the only factor under the control of the operator is the distance 
between the end plates; this is the geometrical restreint the speoimen must 
conform, Therefore, in order to determine the behaviour of the specimen the 
compression stress will have to be plotted as function of the end shortening. 


The unit end shortening, £ , i. 


» the total shortening in спе wave length 
of the shell in axial direction divided by the wave length, oan be easily 
calculated from Eq. (23). It is found that 


2 
Eo E fe ole, d) (s9) 
This equation for the unit end shortening contains only quantities already 
found suoh as the values of f° and + + In Figs. 4 and 5, the compression 
stress & is plotted against the unit end shortening E forg=/ ша 
/120.5 » respectively, It is immediately clear fron this two figures that 
1f the buckling process follows the curves dram, after the shell starts to 
buckle the end shortening has to decrease, In other words, the end plates 
of the testing machine have to move apart, Therefore, the process of 
buokling in this region is highly unstable, since before the operator has 
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times to separate the end plates, the shell will jump to the point P 
(Figs. 4 or 5) which has the same end shortening as the starting point of 
‘the buckling process, but has а much lower compression stress. This jump 
in equilibrium positions involves a release of slastic energy and thus 
explains the rapidity of the buckling process observed and the accompanied 
vibration, 

However, it is still difficult to say whether the shell will jump to 
* state with almost square waves, as shown in Pig. 4, or to а state with 
narrow waves, as shom in Fig. 6. Superficially, one might conclude that 
the narrow waves are the more probable shape because it gives a moh lower 
compression stress, as сап be seen by comparing Figs. 4 and 5. But, 
according to the fundamental concept of mechanics, the true oriterion for 
the most probable equilibrium position compatible with a given value of the 
end shortening is that the elastic energy stored in the shell should be the 
lowest. Ап approximate calculation of the elastic energy of the shells 
shows that for values of )<4/2/, the occurrence of the narrow waves is 
improbable, because the elastic energy stored is higher than that for the 
square waves at same value of the end shortening. However, for “ near 
0,6, and 720/2/ the elastio energy stored in the shell for the narrow waves 
is comparable to that for the square waves at the same value of the end 
shortening. This indicates the possibility of the appearance of narrow 
waves during the very initial stages of the buckling process. 

In any oase, it is certain that there are equilibrium positions of a 
buckled cylindrical shell which involve much lower average compression 


stress Е than that at the beginning of buckling. For instance, in the 


сазе of square wave pattern, the lowest compression stress is given by 


Z= 0/94 £ (40) 
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Incidentelly, this value corresponds olosely to most of the experimental 
results obtained by L. H. Donnell (Ref. 3) and E. E. Lundquist (Ref. 7). 

The corresponding value of the parameter 7 which determines the 
number of waves is equal to 0.225. In case of f 7 1000, the number 
of waves, U, will be 15 which also agrees well with the experimental 
evidence, For this particular value of the radius to thickness ratio, the 
number of waves along the circumference decreases from the 7 = 26 at the 
beginning of the buckling to 7 = 15 at the calculated minimum buckling 
stress. This gradual increase in the sire of waves with the unit shortening 
is also observed by the experiments reported in an earlier paper (Ref. 2). 

Tt is particularly interesting, however, to trace the gradual change 
in the wave pattern during the buckling process. Figs. 6 and 7 show the 
Lines of equal deflection of the wave surfaces corresponding to different 
equilibrium states for two values of the aspect ratio of the wave pattern, 
10 ana A,, These particular equilibrium states are denoted in 
Pigs. 2, 3, & and 5 by а small circle in order to indicate their relative 
position during the buckling process. It is seen that there is а rapid 


shift from the rectangular waves bounded by lines, Ж“ const., and Ра - 


const., as predicted by the classical theory for infinitesimal wave ampli- 
tudes, to staggered rows of circular or elliptical waves, Whereas, the 
rectangular waves are directed alternatively inward and outward, the circular 
or elliptical waves аге ell directed inwards, Тһе transition is practically 
completed for F=4/ Or f, i.e., when the wave amplitude is only 4 ог 6 times 
the thickness of the shell. The occurrence of such inwardly directed 
ciroular and elliptical waves at this stage of the buckling process is in 
good agreement with the experimental observations (Ref. 2). If the experi- 


ment ( Ғ./)0)/ів continued to larger come’ these staggered waves 
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obtain characteristic diamond shapes, The present approximate theory fails 
to give these sharp diamond shaped waves. It is obviously not sufficiently 
exact for such large deflections. Furthermore, when these diamond shaped 
waves occur, the load on the specimen actually falls to в very low value 


such as Ж» 0.06, whereas the theory shows а slight inorease of the 


stress at least for the case ж = 1.0. Therefore, the present oaloulation 
сеп be only considered as а good approximation to the earlier stages of 
buckling when the wave amplitude is only а few times the thickness of the 
shell. Nevertheless, it reproduces the characteristic features of the 
buckling process observed in the laboratory. 


The Effect of the Elastic Characteristic of 


the Testing Machine on the Buckling Phenomenon 


It was stated in the previous paragraph that the state of the specimen 
is determined by the distance between the end plate and that this distance 
is the independent parameter controlled by the experimenter. This statement 
Ав correct only insofar the elasticity in the mechanism of the testing 
machine is neglected. Since there is always a certain amount of elastic 
deflection in the loading mechanism and this deflection is а funotion of 
the load. Hence, if, for example, the loading crank is held at a certain 
position, the compression force acting on the specimen will force the end 
plates apart and thus reduce the amount of end shortening of the speoinen, 
The actual shortening is determined by the load-deflection characteristios 
of the specimen and the testing machine. Assuming that the testing machine 
has а linear elastic characteristic, the compression load is related to the 
end shortening by parallel straight lines, each line corresponding tos зау, 
a constant number of turns of the loading crank. If the loading crank of 
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the machine is held at & fixed position, corresponding values of the 
compression load and end shortening of the speoimen must lie on the straight 


line for this orank position, If the load-end shortening characteristics 


of the specimen itself are given, it is evident that the equilibrium positions 


of the entire system are determined by intersections of the curves re- 
presenting the characteristics of the specimen with the straight lines 
representing the characteristios of the machine, 
Fig. 8 shows representative curves for the characteristics of the 
specimen and two families of straight lines representing the characteristics 


of two différent testing machines. It is evident that after the maximum or 
initial buckling load is reached, the shell will jump to а new equilibrium 
position involving much lower compression load. But this new equilibrium 
position is determined not only by the load-end shortening relationship and 
also by the elastic characteristio of the testing machine. А more elastio 


machine will give & 


t of characteristic straight lines with smaller slopes. ( 
Therefore, in case of curve А (Fig. 8), а more elastio machine will make 

the shell to jump to а higher load, while in oase of curve B, в more elastic 

machine will make the shell to jump to & lower load. This influenoe of the 

elasticity of the testing machine has been discussed by the senior author 

in connection with the plastic buckling of colums (Ref. 8). 


4 Concluding Remarks 
1n the previous paragraphs, the authors bave shown that there are 
equilibrium positions in the buckled shape involving much lower load than 
the buckling load predicted by the classical theory, and thus if the 
specimen is slightly imperfect, it is reasonable to expect much lower 


buckling loads. They have also pointed out that the elastic characteristic 
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of the testing machine might have quite a large influence on the buckling 
process and this might be another cause of the large scattering of the 
data obtained by different experimenters, However, due to the complexity 
of the problem, the results given in this paper can be only considered as 
^ rough approximation and most of the statements made are qualitative 
rather than quantitative, To put the new theory on а 80114 footing, а more 
Acourate solution of the differential equations of equilibrium is necessary. 
Particular attention must be given to the calculation of the elastic energy 
stored in the shell, because it is found that the most probable equilibrium 
depends on the magnitude of the elastio energy stored in the various 
equilibrium positione compatible with the constraint exerted by the loading 
process. 

Furthermore, an inquisitive mind will, perhaps, be pleased by a 
rigorous proof of the validity of all the large deflection equations, 
These equations are established by intuitive arguments, not by systematic 
reasoning. For instance, due to the appearance of sharp curvatures in the 


diamond shaped wave surfaces at large deflections, it is not certain whether 


the curvature of the shell has to be calculated more accurately by taking 
into acoount the second order terms, or the extensions of the median surface 
should be more accurately determined, It is the belief of the authors that 
an investigation of these problems by starting from the general non-linear 
theory of elasticity developed by G. Kirohhoff, J. BoussiWbsq and others 

is very desirable, The recent work by R. Kappus (Ref. 9) is a noteworthy 
contribution in this field of investigation, The senior author has already 
expressed this opinion in his 1939 Gibbs Lecture (Ref. 10) given before the 
American Mathematical Sooiety. 


Section 6 


Take-off from Satellite Orbit 
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Take-Off From Satellite Orbit 
Н. S. Tsient 


Jed ed enne f de AAA с.қ, dosi. АН {жу 


The mass ratio or the characteristic velocity for the take-off of 

a space ship from the satellite orbit is computed for two cases: the 
radial thrust and the circumferential thrust, The circumferential 
thrust is much more efficient in that the required mass ratio is much 
less than for the radial thrust. Both cases show, however, an increase 
of the required mass ratio and the characteristic velocity with a 
reduction in acceleration. With circumferential thrust, the character- 
istic velocity increases by a factor of two when the acceleration is 
reduced from 1/2 g to 1/3000 g. 


συ ΟΙ ΛΠ 


тіс 


For take-off of а rocket from the earth surface, it is convenient to 
have the initial trajectory in the vertical direction and then the thrust should 
e meres |) қалі кә # pue as nat den 

be considerably Art: thas αἵ the rocks Depending upon. 
the relative magnitudes of the aerodynamic drag and the weight, the initial 
ratio of the thrust and the weight should be between 2 and 3 for minimum 
expenditure of the propellant. The situation is quite different for a space 
ship taking off from the satellite orbit: In a satellite orbit, the gravitational 
attraction is completely balanced by the centrifugal force, and the vehicle is 
effectively in a weightless state. This fact has led many fanciers of inter- 
planetary travel to conclude that take-off from satellite orbit requires only 
a very minute thrust. For instance, L. Spitzer (Ref. 1) proposed a nuclear 
power plant for з space ship to be accelerated at only 1/3000 g. Another 


example is the extensive discussion of interorbital transport techniques by 
н. Preston ~ Thomas (Ref. 2) based upon the assumption of equally small 


acceleration. On the other hand, W. von Braun (Ref. 3) seems to prefer а 
very much larger acceleration of approximately 1/2 g for take-off from the 
satellite orbit. 

The magnitude of the acceleration has a strong bearing on the optimum 
type of power plantto be used: The ion-beam rocket is only possible for 
very small acceleration, while for moderate acceleration, chemical rocket 
is required. Therefore the question of the magnitude of acceleration is an 
important one for interplanetary flight. The purpose of this note is to 
x relation between the acceleration and the mass ratio 
required for escape from the earth's gravitational field,starting from the 
satellite orbit. It is hoped that the present investigation will give the future 
generation of astronautical engineers a mere rational basis for designing 
space ships. 


Вавіс tions 

The problem considered is the motion of a space ship under the influence 
of the rocket thrust and the gravitational attraction of a single massive body, 
say the earth. Then if the rocket thrust is in the plane of trajectory, the 
trajectory of the space ship will remain ina plane. Let the position of the 
ship at any time instant t be given by the polar coordinates г and ў : 三 
the distance from the center of attraction, and ϐ the angular position. If 
the components of the rocket thrust per unit mass of the vehicle are В in the 
radial direction and Ө in the circumferential direction, and if g is the 
magnitude of gravitational attraction at the starting satellite orbit = 7; 


(Fig. 1), then the equations of motion of the space ship are 


45- Rs "Ш-Н 


ЖҮ #)- ro a 
By using the subscript 0 to indicate quantities at the starting instant 
%-0, the equilibrium condition of the satellite orbit is given by 


ШЫР o 


Initially, the radial velocity is zero, i.e., 


[77] =0 (4) 


These аге the initial conditions. Жоғ the space ship to have sufficient energy 
to escape the earth gravitational field at the end of the powered flight, the 
sum of the kinetic energy and potential energy must vanish at the end of the 
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end of the accelerating period. Let that instant be denoted by the subscript 1. 


2/8) 4 2 

ҮЛ * FE) 7 ZË E 6) 
With апу specified variation οἵ the thrust forces R and @ ав functions 

of time, the above system of equations determine completely the take-off 

trajectory of the space ship. In the following sections, two special cases of 

practical significance will be discussed in detail: the case R = constant, 


Thus at t=t, 


© #0, purely radial thrust; and the case R= O. O = constant, purely 
circumferential thrust. 


Radial Thrust 
If the thrust is always radial and is propertianal to the instantaneous 
mass of the vehicle, а nonfimensional thrust factor #1 can be introduced as 
5 c 
Furthermore, let 
e 7 


š ةت اا م‎ жаі al 
Eqs. (1) ава @ аъ 


е ΠΝ Á 
£(r#)=o е 


4 
Eq. (8) сап be immediately integrated and by using the initial condition of 


Eq. (3), the result of integration is 
m" 


4 
LE 
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d 
By substituting this equation into Eq. (7), the final equation for ¢ ів 


-rH 


The normimensiona1 radial velocity is df/dt. This is related to the 
physical radial velocity d r/dt as follows 


PE du 
Eq. H9) can be rewritten as 
HI nan 

Since d(/dt = 0 when 1 = and gat according to Eq. (4), the result of 
integrating the above equation is 

TS 

fit) m 4p 18-1) hi d) ite 2) 
Therefore the non-dimensional time Ф can be calculated as a function of the 
radius f as follows 


ПЕ Rm s n 
ШАРЛЫ f£) 
With Eqs. (9) and (Ê), the end condition of Eq. (5) can be written as 


ΑΛ ΗΡΙ Hege 


Or simply 


و 


Then the velocities at the end of acceleration period are 


7 


ГАР) = ау Hr age 


- E 
The time у for the powered flight can be obtained from Eq. (15) by setting 
the upper limit of integration to $}. * f Ας ав * 


MEE ыы ыы 


where F and E are the elliptical integrals of first kind and second kind 
respectively. 

If M(t) is the instantaneous mass of the space ship, and с the effective 
exhaust velocity of the rocket, then 


tgs cd . 


Therefore the mass ratio м/м can be calculated as follows: 


Нм) = WEE per, 


By using the result of Eq. (16), 


ΜΑ AE ) 14 
die ως МИЎ. used аже а j өт 


ra 
M 
When the acceleration is very large, />> 1, the integrand in Eq. (15) can 
be expanded in terms of this parameter. Then the mass ratio is (М, αἱ 
€ 
‚= а nee s Ms 
pr el) Ша v Ж ae) 
Y 
The relation of Eqs. алам is plotted in Fig. 2. For е ж 4 
the mass ratio becomes infinite, The reason is that at this value ой 


* The author is indebted to Dr. Y. T. Wu who kindly supplied the relation 
of Eq. ep. 
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acceleration, there is a radial position where the thrust force is equal to 

the gravitational attraction and ome further increase in the energy of the 
vehicle can occur. Therefore the radial thrust per unit mass, if maintained 
constant throughout the powered flight, should be larger than 1/8 g. With 
increasing thrust, the required mass ratio for escape from the earth's 
gravitational field decreases. VThe asymptotic value of Ань fu ) is ipe 1 
However there is no appreciable improvement in going to higher thrust than 
18. This strong ος of the mass ratio upon the acceleration factor" 


ent is contrary to opinion that for take-off from satellite orbit only very small 
thrust is required. 一 mMm 


1f the thrust is always circumferential and proportional to the mass of 
the vehicle, then а nas Hust e je Ç tan dt CK enk Kab 


©» vp 
By using the same non-dimensional variables as defined in Eq. (6), the 
equations of motion are 


Ae s) w 


‘The initial conditions of Eqs. (3) and (4) are 


Ші, ize, € get, to 


Therefore Eq. eu gives another initial condition that 


HAE) =° 


] — 
By eliminating Í from Eqs. (29) and (23), 


f dk ee v, 

This is a third order differential equation with three initial conditions 
specified by Eqs. (23) and (21). No simple general solution can, however, 
be obtained. 0000 а 
that are valid for large values of ¢ от рожи values of y . 

For very large values of y , the acceleration period is expected to 
be short and the change of the radial position to be small. Then the value of f 
must be very close to the initial value of unity, By taking p to be unity, 


Eq. (24) becomes εί... : 


7 aC rr + уе 


беге C is the Negration 3 C, however, must be 1 because of 
the initial condition of Eq. (2h. The appropriate approximate solution for 
4 for very large V is thus 
б) 37 
* /+ ўт? + уч z" . өз 
To obtain higher terms in this power series, the usual series substitution 
method may be used. The calculation is somewhat lengthy and therefore 
will not be reproduced here, The result is 
E 


= 1+ әт» dvr Zr Bes! 4--- (26) 


By using the result of Eq. αὖ, the radial velocity is obtained by 
differentiation. Then Eq. (28) gives the circumferential velocity. The end 
* condition of Eq. (5) can be modified into the following more convenient form 
by multiplying it yar 


o= [б Na 
By substituting the solution of Eq. (26) into this condition, an equation for 
determining 7, is obtained: 


2 ov 
0= / Mte ty EY УТ, t glnA Els) 4. n 


The mass ratio /,/M, can be calculated in the same way as in the 
previous section and can be determined through the new parameter x defined as 
follows: 


ipa ^g Inm ή-α 


Eq. eh then can be written as 


ё 


0= M= r r n 


өй у S ا‎ 
ТМ жүз 
Since the calculation is designed for large values of Y , the appropriate 
expansion of x'should be a series in inverse powers of У. Eq. (88) 
suggests specifically 


x? 
у) = x". 4 ipee 


V. e. 
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where KD XU) naa xU) aze constants independent of V . By 
substituting Eq. (38) into Eq. (89) and equating equal powers of y , the 
following set of equations results: 
* 
* 09-7 0 бт) 


42 2 
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The explicit numerical solutions are then 


"= 2 =1 = Ma] * 
а en 
X = 6.062349 

* Ace 


This completes the calculation of mass ratio for large values of the 
acceleration factor у. 
For the other extreme case of very small values of f , it is to be 
expected that fn Eq. (АУ) the acceleration will be very smallandithe rm 
rh, will be very much smallbr than f . Therefore a good approximation 
( 
of Eq. (24) at small у is 


w Ix 4 


The solution of this equation with the initial condition of fiat 7-04 


9= on 


4 * 37 
LL E: 
At 7-0 , the radial velocity and the radial acceleration are thus not zero 
as required by the initial conditions of Eqs. ай and αὖ. They are however 
very small, because У is very small. Therefore the solution of Eq. eh 
is a good approximation to the exact solution. 
To the same approximation, Eq. (20) becomes 
34 


rž- ж-< f1-vz) (m 


-- ^ mall = 
( This means thatipecause of the extremely slaw acceleration, the сетті му 
\ force per unit mass T y practically balances the gravitational 
attraction gt every instanD The end condition of Eq. (5) can then be written 


as 


ye pa 54 
dee - (а-о 65 


where X is үр . The appropriate solution for X is then 
A 


Х= /一 pif 


Since the mass ratio 22 , is related to X by Eq. αὐ, Еа. dn 
actually gives the mass ratio for escaping the gravitational field with very 
small acceleration. : 
The parameter Ж is plotted against V in Fig. 3, using both Eq. (30) 
with Eq. (34) and Eq. (48). when У approaches sero, X approaches 1. 


ЖҰЖ 
th andiod ponte 


a X. he 
« Ар ue 
y bera 4. At 
b 


pe 


hence the mass ratio / М, , decrease monotonically. Therefore, same as 
the result for purely radial thrust, there is а strong influence of the magnitude 
of acceleration on the required mass ratio. However as far as decreasing 
the mass ratio is concerned, there is no appreciable advantage in using У 
greater than 1/2. 
— я 
Discussion 
By comparing Fig. 2 with Fig. 3, it is apparent that the radial thrust 
is much less efficient than the circumferential thrust for take-off from the 
satellite orbit. For large thrusts, the value of log (/% / М) for radial thrust 
is more than twice that for circumferential thrust. Furthermore, in case 
of radial thrust, the ratio of thrust to the instantaneous mass, if maintained 
constant, must be larger than /e. In case of circumferential thrust, no " 
such limit exist, Therefore circumferential thrust is definitely preferred. 
The quantity ¢ log /M,/M,) is a measure of the performance or the 
capability of the vehicle. It has the dimension of a velocity and is actually 
the increase of velocity which the vehicle is capable of in space without 
gravitation. This quantity is conveniently called the characteristic velocity 
of the vehicle. Let this be denoted by У. Then for the case of 
circumferential thrust, Eq. αὖ gives 


V= fA n) ря я % 


where (| is the “escape velocity" from the surface ofthe earth and А is 
the ratio of the radii of the satellite orbit and the earth, „ is equal to 
+ — 1Mj2 km/sec. Fig. 3 then shows that by decreasing the acceleration 
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1/2 to 1/3000 g, 7 , hence the required characteristic velocity И, will 
increase bya factor of two. This is a very important point for the 
designers of space ships. 
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Circumferential Acceleration 
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